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The concept of self-organized criticality (SOC) 
concerns the spatiotemporal complexity in the sys- 
tems that contain information over a wide range 
of length and time scales [1,2]. It implies that 
through a dynamical process a system can start in 
a state with uncorrelated behavior and end up in 
a complex critical state with a high degree of cor- 
relation. This concept and the prototype model, 
sandpile model, are proposed by Bak, Tang and 
Wiesenfeld [3] in 1987. Self-organized criticality is 
so far the only known general mechanism to gen- 
erate complexity [4], and hence the one trying to 
understand why nature is complex, not simple. 

Systems which can exhibit SOC are common in 
physics, geography, biology, and even social sci- 
ences such as economy. Such kinds of complex 
phenomena are ubiquitous in macroscopic world. 
Recently, it has been proposed by Meng et al. [5] 
that SOC may exist in microscopic systems — at 
the level of quarks and gluons, as well as in macro- 
scopic world. 

Evidence from the biology has suggested [6,7] 
that the ecology of interacting species had self- 
organized to a critical state. In 1990, Bak, Chen, 
and Crutz [8] created a cellular automaton sim- 
ulating a society of living organisms operates at, 
or very close to, the critical state when driven by 
random mutations. However, the model is very 
sensitive in the sense that small modifications of 
the details may drive it away from the critical 
state. NKC model proposed by Kauffman and 
Johnsen [9] can exhibit a transition from order to 
disorder, but the criticality emerged in the system 
is obtained through parameter tuning, not self- 
organizing. In 1993 [10], a simple model of evolu- 
tion, Bak-Sneppen model, was introduced by Bak 
and Sneppen. Instead of considering the evolu- 
tion on the individual level they present the co- 
evolution of species on a coarsed-grained scale. In 
their model the whole species is represented by a 
single fitness, i.e., a random number chosen arbi- 
trarily from a uniform distribution between zero 
and 1. And mutations correspond to updating the 
fitnesses of a given extremal site and its two near- 
est neighbours with three new random numbers 
chosen from the same flat distribution between 
zero and 1. Such model of an evolving biology 
can self-organize to a critical steady state during 
which all sizes of avalanches can occur. Most im- 
portant of all, their model can exhibit punctuated 
equilibrium behavior observed in biology. Two ex- 
act equations describing self-organization and the 
average avalanche size behavior respectively were 
found later [11,12]. A hierarchical structure of 
avalanches were also observed in B-S model. It 



was even proposed by Ref. [2] that the formation 
of fractal structures, the appearance of j noise, 
diffusion with anomalous Hurst exponents, Levy 
flights, and punctuated equilibria can all be re- 
lated to the same underlying avalanche dynamics. 

Our model of evolution is based on B-S model, 
but differs decisively in driving mechanism of in- 
teraction between neighboring species. We also 
consider coevolution of an interacting species sys- 
tem and each species is represented by a single 
fitness, i.e., a random number chosen arbitrarily 
from a flat distribution between zero and 1. But 
when considering the interaction between neigh- 
boring species we impose a constraint concerning 
the differences of fitnesses between that of the ex- 
tremal species and those of its nearest neighbours 
at each time step. Before knowing how the con- 
straint is imposed let us take a first look at another 
case of evolution, a non-interactive biology. In a 
non-interactive biology each species would tend to 
evolve towards a stable state where the fitness of 
each species approaches 1, but the evolution pro- 
cess is extremely slow. If we use a/ to denote 
interaction strength, which represents the degree 
of interaction between the extremal site and its 
nearest neighboring species, it is natural to let cti 
be 1 in B-S model and a/ be zero in the non- 
interactive biology. If so, it seems that these two 
cases of ai correspond to two extremal cases of 
interaction . Then, if ai is allowed to take any 
value between zero and 1, what can we do with 
our model? Thus, several questions arise there: 

1) What does interaction strength mean for an 
evolution model? 

2) How to present the definition of interaction 
strength aj and impose it into a model of evolu- 
tion? 

3) If the two questions are solved, then are the 
model of evolution and some features of it affected 
by the values of interaction strength? 

It will be shown that these questions could be 
solved successfully in the following text. 

Our model is intended to consider the evolu- 
tion of an ecosystem which consists of a number 
of species. Followed the ideas proposed in Ref. 
[10] each species is represented by a single fitness. 
The fitness may represent population of a whole 
species or living capability of the species. A high 
fitness of a certain species may imply that number 
of the species is immense or its living capability 
is very great, and vice versa. Change of fitness 
of a species may imply change of number of the 
species or change of its living capability. So it is 
natural to expect that a species with a low fit- 
ness is more likely to change, namely mutation in 
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biology, in order to live better and/or longer in 
nature. Only through mutations the species with 
a bad fitness can have the chances of choosing a 
better fitness in order to avoid extinction. Fur- 
thermore, the fitness of each species is affected by 
other species which are also parts of the ecosystem. 
Any adaptive change of any species may change 
the fitness and the fitness landscape of its coevo- 
lutionary parts coupled in the same ecosystem. So 
the species may interact with each other through, 
say, a food chain. Hence, the species with a high 
fitness may live well and comfortably unless its 
bad neighbours are going to mutate. 

Our model is defined and simulated through the 
following items: 

(1) A number of, say, L d species are located 
on a d-dimensional lattice of linear size L. Ini- 
tially, random numbers chosen arbitrarily from a 
uniform distribution between zero and 1, p(f), are 
assigned independently to each species. At each 
time step, 

(2) choose the extremal site, that is, the species 
with the lowest fitness, / m i n , among all the species 
and update it by assigning a new random number 
also chosen from p(f) to it and 

(3) mutate those of its 2d neighbours whose fit- 
nesses satisfy the constraint /, — / m i„ < a by as- 
signing new random numbers between zero and 1 
to them (/j denotes the fitness of the ith nearest 
neighbours.), a is a parameter between zero and 
1 and is fixed for a given model. 

It should be pointed out that different values 
of a correspond to different versions of the model, 
even to say, different models. Consider two special 
values of a, and 1. For a — the model returns 
to the non-interactive biology since difference of 
fitness between those of any two neighbours is al- 
ways greater than zero, so in the model in which 
a = none of the neighbours of any extremal 
species will be chosen for updating at each time 
step. This is the case where there is no interac- 
tion between neighbours. For a = 1 the model 
returns to a <i-dimcnsional B-S model. It is be- 
cause that difference of fitness between those of 
any two neighbours is less than I so at each time 
step all the 2d nearest neighbours of the extremal 
site will be chosen for updating. 

For a given a let the model evolve from the be- 
ginning when the first extremal species is chosen 
and we determine how many of its 2d neighbours 
will be chosen, according to the constraint on dif- 
ference of fitness, for updating at the same time 
step. The updating process, i.e., the evoultion of 
the model, continues indefinitely. For a given a 
satisfying < a < 1 the updating of the 2d near- 



est neighbours of the extremal site at each time 
step will have many kinds of probabilities. For 
instance, maybe none of the 2d neighbours will 
be updated, maybe half of the 2d neighbours will 
be chosen for updating, etc. If we do not distin- 
guish the neighbours when we only care the num- 
ber of the updated neighbours at each time step it 
is straightforward that the updating of the neigh- 
bours will have 2d + 1 probabilities. Say, at some 
time step, an extremal site and m of its 2d near- 
est neighbours are chosen for updating according 
to the constraint. Such an updating is called by 
us a m-event. As the evolution goes on we can 
observe many kinds of such events, during which 
m can be different. If the evolution time is large 
enough we can obtain 2d + I kinds of events dur- 
ing which m spans through and 2d. If we define 
Pd{m) the probability of m-event among all events 
during time period T, i.e., 



Pd(m) = 



N(m) 



(1) 



where N(m) is the number of m-event during the 
time period T and Nt, the total number of all 
events, and d denotes dimension. In the T limit 
(T 3> L d ) wc will obtain the distribution of m- 
event, that is, 



Pd{m) = lim 



N(m) 
N T ■ 



(2) 



And Pd(m) should satisfy the normalization, 



=2(2 



(3) 



m=0 



Next we will present the definition of interaction 
strength ai related to P<j(m), 
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2d 



^Yd^ mPd{m) - 

m=0 



(4) 



One can see that ai is actually the statistical ratio 
of number of updated neighbours among 2d ones 
during the evolution. This can be easily seen from 
two extremal cases when a = and a = 1. For 
a = 0, ai = 0, this is because Pd(0) = 1 and 
Pd{m) = for < m < 2d; for a = 1, a/ = 1, 
this is because Pd(2d) = 1 and Pd(m) = for 
< m < 2c?. That is, in the nointeractive bi- 
ology interaction strength ai is zero while in the 
B-S model ai is I. Hence, it is natural to ex- 
pect that the definition of ai can give a good 
description of the change of strength of interac- 
tion between neighboring sites when a changes. 
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So, it is also natural to expect that < aj < 1 
for < a < 1, and that a/ increases as a in- 
creases. There should also exists one-to-one cor- 
respondence between aj and a. Having these in 
mind we measure the distribution Pd{m) and ai 
for ten different values of a for one- and two- 
dimensional models on the computer. Simulation 
results are given in Fig. 1 and Fig. 2 respectively. 
Figures (a) and (b) in Fig. 1 show the distribu- 
tion of m-event (0 < m < 2d) for ten different 
values of a in one- and two-dimensional models 
respectively. It is clearly shown in the above fig- 
ures that Pd(0) decreases with increase of a and 
P d (2d), i.e., P d=1 (2) for d = 1 and P d=2 (4) for 
d = 2, increases as a increases. Note in the fig- 
ures change of Pd{m) (0 < m < 2d) with increase 
of a. Their exists a peak in the curve of a depen- 
dence of Pd(m) (0 < m < 2d). This can ensure 
the normalization of Pd(m). Two plots in Fig. 2 
show the dependence of ai on a in one- and two- 
dimensional models respectively. As shown in Fig. 
2 ai almost increases linearly as a increases, and 
most important of all is that one aj corresponds 
to one and only one a. From here one can clearly 
see that our definition of ui can be explicitly re- 
lated to a and hence can be put into our model 
naturally. Relating Fig. 1 and Fig. 2 one can 
see that increase of interaction strength increases 
P d (2d). That is, when ai increases any extremal 
site will possibly affect its nearest neighbours more 
strongly. Furthermore, different values of ai cor- 
respond to different interaction strength. Gener- 
ally speaking, larger ai represents stronger inter- 
action, and vice versa. Hence, the definition of ai 
can provide a good description for the strength of 
interaction between neighboring sites. It is a good 
quantity in describing the interaction. In the fol- 
lowing text dependence of self-organized threshold 
and some critical exponents on cti will be given. 
Since ai has one-to-one correspondence with a 
and increases as a increases, it is hence convenient 
and equivalent to present the dependence of these 
quantities on a. 

The model is already defined, it is natural to in- 
vestigate whether the model can self-organize to a 
critical state. That is, we should observe the "fin- 
gerprint" of SOC [4]. If so, it is also worthwhile 
to know whether the criticality is sensitive to the 
value of a. In addition, the self-organization, re- 
ferred to a dynamical process whereby a system 
starts in a state with uncorrelated behavior and 
ends up in a complex state with a high degree of 
correlation, of the system and punctuated equi- 
librium, the most important feature of evolution, 
should also be observed. 



Fig. 3 shows the space-time fractal activity pat- 
tern for a one-dimensional model of size L = 100 
with a — 0.8. We track the updated sites at each 
time step. S and R in the figure denote num- 
ber of updated time steps and location of updated 
sites respectively. Simulations of one- and two- 
dimensional models for different values of a are 
also done and exhibit the similar space-time frac- 
tal activity. Indeed, spatiotcmporal complexity 
emerges in our model and its appearance is in- 
dependent of the value of a choosen. 

We now present some explanations of the quan- 
tities which will appear in the following equations 
for those readers who are not so familar with SOC. 
/min denotes the extremal fitness at each time step 
in the evolution. G(s), the gap appeared in punc- 
tuated equilibrium, is an envelope function that 
tracks the increasing peaks in f m [ n . Its defini- 
tion is : at time step s the gap G(s) is the maxi- 
mum of all the minimum random numbers chosen, 
/min(s'), for all < s' < s. f c is the value of G(s) 
at critical state, i.e., 

f c = lim G(s). (5) 

Sg(s) is the size of avalanches correspond to pla- 
teus in G(s) during which ,/ m i n (s) < G(s) and 
(S)g(s) is the average value of Sg( s ) ■ An avalanche 
is denned as subsequent mutations below a certain 
threshold. Hence with this definition there is a 
hierarchy of avalanches each defined by their re- 
spective thresholds. So we can have /o-avalanche 
where /o is only an auxiliary parameter between 
zero and 1 to define avalanches. More detailed 
definition of /o-avalanche is given in Rcf. [5] . The 
size of an avlanche, S, is the number of subse- 
quent mutations below the threshold. And 7 is a 
critical exponent which governs the divergence of 
Sg(s) when s approaches infinity. And 

^cov IS the 

number of sites covered by an avalanche. Appar- 
ently n cov < (2d+ 1)S (S is the avalanche size) in 
d-dimensional space. (n cov ) denotes the average 
value of n cov . The above defined quantities will 
have their counterparts in our model. So we will 
use the same definitions of these quantities while 
make some minor corrections on the symbols of 
them. 

Following the method used in Rcf. [2] we moni- 
tor the extremal signal f mm as a function of s dur- 
ing the transient in the one- and two-dimensional 
models for different values of a. Again, we ob- 
serve Devil's staircase [4] in all these cases. Fig. 
4 shows punctuated equilibrium behavior in one- 
dimensional model of size L = 100 with a = 0.5. 
Hence, we can see punctuated equilibrium does 
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emerge in our model of evolution and its emer- 
gence is a independent. 

Observations through above simulations suggest 
that our model can self-organize to a critical state. 
But how to determine the self-organized thresh- 
old f c is still a hard bone to us. Fortunately 
Ref. [2] provides us a method. In Ref. [2] 
the author presents two exact equations for Bak- 
Sneppen model. With a reasonable scaling ansatz 
of average size of avalanches they can determine 
f c very accurately. It is straightforward to ex- 
pect that the two exact equations found in Bak- 
Sneppen model have the corresponding ones in our 
model. According to the derivation of these two 
equations [2] we can directly write down two sim- 
ilar exact equations in our model except for some 
minor replacements on the symbols of some quan- 
tities. Introduction of interaction strength ai sug- 
gests that some quantities are ai and hence a de- 
pendent, which can be seen from our simulations: 
as a increases f c has a tendcy to decrease. Having 
this in mind we make replacements of symbols of 
some quantities: 



G(s) 

fc 

7 
r 
D 



G(s,a) 
■* fc(a), 
-> 7(a), 
-> r(a), 
-> D(a). 



Two exact equations are given as below, 
dG(s,a) _ 1 - G(s,a) 



ds 



L d (S) 



G(s,a) 



dln (S)/o _ ( 



dfo 



l-/o 



(6) 



(7) 



In order to solve the two equations a scaling 
ansatz of (S)Q( SyCt } should be given, 

<S}G(s, Q )~[/e(a)-G( S ,a)r 7(a) . (8) 
Inserting Eq. (8) into Eq. (7) one obtains 
n cov ) f [fc{a) - G(s,a)] 



7(a) = lim 



1-/0 



(9) 



Using Eq. (g) we can determine f c (ce) and 7(a) 
very accurately. We measure f c {a) and 7(a) in 
the one- and two-dimensional models for different 
values of a. It is to our expectation that these 
two quantities are a dependent. Fig. 6 shows the 
dependence of f c { a ) on a an d figure (a) in Fig. 7 
shows that of 7(a) on a. In addition, we measure 



two basic exponents of the model, T(a), avalanche 
size distribution, and D(a), avalanche dimension. 
We find they are also a dependent. Figures (b) 
and (c) in Fig. 7 show their dependence on a 
respectively. 

The following text will present our analysis 
of the simulations and our conclusion. Let us 
first check the "fingerprint" of SOC. Fig. 5 
shows avalanche distribution P ave d(S) for one- 
dimensional model of size L = 100 with a = 0.7. 
Such distribution can also be found in one- and 
two-dimensional models for different a. Indeed, 
power law emerges in our model and its appear- 
ance is a dependent. Then let us come to the 
results of f c {ot)- Two plots, (a) and (b), in Fig. 6 
show the dependence of f c (ct) on a in the one- and 
two-dimensional models respectively. It is clearly 
shown that these two figures display the similar 
behaviors of f c (a). Firstly, we can see that in 
both figures f c (a) decreases as a increases. It is 
not difficult to understand this kind of behavior. 
As a increases, i.e., ai increases, the chances for 
the nearest neighbours of a given extremal site to 
be chosen for updating at each time step will be 
greater. This can be easily seen from Fig. 1. As 
shown in Fig. 1 Pd(2d) increases as a increases 
and reaches 1 for a = 1, and Pd(0) decreases as a 
increases and reaches zero when a = 1. So, with 
the increase of a , that is the increase of ai , more 
possible neighbours will be involved in the evolu- 
tion, hence, the threshold / c (a) will be lowered 
further. This can be explained in another way. 
Compare the values of f c (a) for the same a in 
one- and two-dimensional models. Say, compare 
/c(l) in one dimensional model with that in two- 
dimensional model. It has been shown that / c (l) 
in the former case is greater than that in the latter 
case. In one -dimensional model in which a = 1 
when an extremal site is chosen, two of its near- 
est neighbours will also be chosen, for updating at 
each time step. While in two-dimensional model 
with a = 1, an extremal site, together with its four 
nearest neighbours, will be chosen for updating at 
each time step. Why we mention this two cases 
here is just trying to show that the increase of av- 
erage number of neighbours involved in the evolu- 
tion will lower the self-organized threshold. Thus, 
one can see, the increase of interaction strength 
will increase the number of neighbouring sites in- 
volved in the evolution and hence lowers down the 
value of f c (oi). Secondly, it is shown in both figures 
of Fig. 6 that / c (a) almost decreases linearly as 
a increases for < a < 0.6 and decreases asymp- 
totically as a increases for a between 0.6 and 1.0. 
This implies that effect of the constraint on the 
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evolution grows implicitly as a increases. When 
interaction strength is small the effect is very ex- 
plcitly. But as interaction strength increases the 
effect will be not so explicitly shown. Specifically 
we measure f c {ct) when a = 1 for one- and two- 
dimensional models. We find / c (l) = 0.668±0.001 
for d = 1 system of size L = 100 , which is very 
close to the value in Refs. [2] and [13], who found 
/ c (l) = 0.660702 ± 0.00003; for d = 2 system of 
size L = 20, / c (l) = 0.334 ± 0.00006 which is 
close to the corresponding value in Rcf. [2], who 
found / c (l) = 0.328855 ± 0.000004. In addition, 
for d = 1 we find / c (0.1) = 0.96388 ± 0.000002 
for a system of size L = 100 and for d = 2, 
/c(0.1) = 0.9179 ± 0.00001 for a system of size 
L = 20. 

We also measure 7(a) for different a in one- 
and two-dimensional models. It is found that 7 
is also a dependent. Dependence of 7(a) on a is 
given in Fig. 7. In figure (a) in Fig. 7 7(a) first 
increases and then decreases as a increases. For 
d = 1, 7(1) = 2.6166 ± 0.00004, which is close 
to the value found in Refs. [2,13,14], who found 
7(1) = 2.70 ± 0.01. For d = 2, we find 7(1) = 
2.249 ± 0.0004, which is not in agreement with 
Rcf. [2], who found 7(1) = 1.70 ± 0.01. Maybe 
this is because of our small size of system. If the 
size is larger the results will be more precise. 

For different values of a in one- and two- 
dimensional models we measure two basic expo- 
nents, r(a), which characterizes the distribution 
of avalanche sizes, and D(a), the avalanche dimen- 
sion. The definition of r(a) is P ava i(5) - S- T< - a \ 
and that of D(a) : n cov ~ S D< - a ^ d . Results of 
these two basic exponents are given in Fig. 7(b) 
and Fig. 7(c) respectively. It is shown in the plots 
that r(a) first decreases as a increases and then 
changes slowly with variation of a. Specifically 
for d = 1, we find r(l) = 0.858 ± 0.001, which 
is not in agreement with Ref. [2], who measured 
t = 1.07 ± 0.001, but agrees with Rcf. [10], who 
measured r = 0.8 ± 0.1. This is because our size 
L = 100 is close to the size L = 64 in Rcf. [10] but 
far from the size L = 10 4 in Rcf. [2]. For d = 2, 
r(l) = 1.131 ± 0.0005, which is close to Ref. [2], 
who measured r = 1.245 ± 0.01. 

Fig. 7(c) shows the dependence of D(ot) on a 
for one- and two-dimensional models. It can be 
inferred from the plots that D(ot) first decreases 
rapidly and then slowly as a increases. For d = 1, 
D{1) = 2.4189 ± 0.0001, which is close to Ref. 
[2], who measured D = 2.43 ± 0.01, and for d = 
2, D(l) = 2.94 ± 0.0001, which agrees with Refs. 
[2,9,10], who measured D = 2.92 ± 0.02. It is 
clearly shown that for all a D(ot) is larger than 



the dimension space d. 

It should be emphasized that interaction 
strength, oti, is not a parameter tuning the model 
to a critical state, since in our model when a model 
is chosen the corresponding ai is fixed during the 
whole evolution process. Thus, criticlity emerged 
in our model is self-organized, not tuned. And ap- 
pearance of criticality in our model is independent 
of ui chosen, which you can test on your own PC. 
This can strongly support the idea that SOC does 
not depend on the dynamical details of the system. 
That is, self-organized criticality should be univer- 
sal. Furthermore, self-organized fractal growth is 
basically different from growth processes, say, de- 
scribed by (variants of) the Kardar-Parisi-Zhang 
(KPZ) equation [15,16]. The KPZ equation is 
scale invariant by symmetry, thus the criticality 
is not self-organized. SOC cannot, even in princi- 
ple, be regarded as sweeping a system through a 
critical point, which contrasts to the claims in Ref. 
[17]. SOC should be an attractor of the complex 
system, but this attractor is vastly different from 
the one found in chaos — " strange attractor" . 

As shown in this paper and in others [2,4,10], in- 
teraction plays a very important role in the mod- 
els which exhibit SOC. If there is no interaction 
between individuals in a system, the system will 
evolve towards a frozen state and the evolution 
process will be indefinitely long. This is clearly 
shown in non-interactive biology in which fitness 
of each species tends to be 1. Our model and 
simulations of its different versions corresponding 
to different degrees of interaction between neigh- 
bours imply that only the cocvolutionary system 
can evolve to a self-organized critical state despite 
the fact that interaction strength may be rela- 
tively small. Through our simulations that we 
have learned it would take longer time for the sys- 
tem to involve to a critical state when the interac- 
tion strength oti is smaller. It is also worthwhile 
to perform two simulations of the models in which 
oti is very close to zero and 1 respectively. We can 
expect in the former case the model will evolve 
to a frozen state, while in the latter case we will 
approach B-S model. 

Another important feature of our model is that 
fitness itself is directly involved in the interaction, 
which is not realized in the B-S model. As a co- 
evolutionary system, interaction between any ex- 
tremal species and its neighbours should and must 
involve the features of the extremal site and that of 
the neighbours. Because of model's simplicity each 
species has only one feature: fitness, represented 
by a random number chosen arbitrarily from a flat 
distribution between zero and 1. Hence, this fea- 
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ture should enter the evolution model. Through 
a constraint which can be related to interaction 
strength the fitness is involved in and injected into 
the evolution process of the system. It is shown 
in our model that evolution is indeed an coevolu- 
tionary phenomena, which agrees with Darwin's 
opinion on evolution of the biology [18]. 
Thus, in conclusion: 

(1) A simple model of evolution with interac- 
tion strength defined and considered is proposed. 
The models with different interaction strength 
(0 < ctj < 1) can self-organize to critical states. 

(2) Simulations of one- and two-dimensional 
model of various degrees of interaction strength 
show that f c (aj) decreases as ai increases. It is 
also shown that 7(0;/), and two basic exponents , 
r(ai) and D(ai), are ai dependent. 
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Figure Captions: 

Fig. 1: Distribution of m-event for different a 
in the evolution of (a)one-dimensional models of 
size L = 100 and (b)two-dimensional models of 
size L = 20. 

Fig. 2: Dependence of interaction strength aj 
on a in the evolution of one-dimensional models 
of size L = 100 and two- dimensional models of 
size L — 20. 

Fig. 3: Space-time fractal activity pattern for 
one-dimensional evolution model with a = 0.8. 
S is the number of updated steps and R is the 
location of updated sites. 

Fig. 4: Punctuated equilibrium behavior 
emerges in one-dimensional evolution models of 
size L = 100 with a = 0.5. G(s) is the gap that 
tracks the peaks of extremal signal, / m ; n , in the 
transient. f c in the plot is about 0.8556. 

Fig. 5: Avalanche distribution for one- 
dimensional model of size L = 100 with a = 
0.7. S, the size of avalanche, is the number of 
subsequent mutations below the threshold 0.728. 
-f'avai('S') denotes the distribution of S. The slope 
of the curve is about 0.774 ± 0.001. 

Fig. 6: Dependence of f c (a) on a for one- 
dimensional evolution models of size L = 100 and 
two-dimensional models of size L = 20. 

Fig. 7: (a) Dependence of 7(a) on a for one- 
dimensional evolution models of size L = 100 and 
two-dimensional models of size L = 20. (b) De- 
pendence of r(a) on a for one-dimensional evolu- 
tion models of size L = 100 and two-dimensional 
models of size L = 20. (c) Dependence of D(a) 
on a for one-dimensional evolution models of size 
L = 100 and two-dimensional models of size L = 
20. 
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